The dynamical Casimir effect in a superconducting coplanar waveguide 
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We investigate the dynamical Casimir effect in a coplanar waveguide (CPW) terminated by a 
superconducting quantum interference device (SQUID). Changing the magnetic flux through the 
SQUID parametrically modulates the boundary condition of the CPW, and thereby, its effective 
length. Effective boundary velocities comparable to the speed of light in the CPW result in broad- 
band photon generation which is identical to the one calculated in the dynamical Casimir effect for 
a single oscillating mirror. We estimate the power of the radiation for realistic parameters and show 
that it is experimentally feasible to directly detect this nonclassical broadband radiation. 

PACS numbers: 85.25.Cp, 42.50.Lc, 84.40.Az 



o 
o 
I 

Oh: 



O 

o 



> 

m 

o 

o\ 
o 



X 



Two parallel mirrors in empty space are attracted to 
each other due to the vacuum fluctuations of the elec- 
tromagnetic field, because of the different mode density 
inside compared to outside of the mirrors. This striking 
effect of quantum electrodynamics (QED) was predicted 
by Casimir in 1948 and since then it has also been verified 
experimentally (see, e.g., Ref. [H). 

If the mirrors move, there is also a mismatch between 
vacuum modes at different instances of time. It was pre- 
dicted that this may result in the creation of real pho- 
tons out of vacuum fluctuations. This dynamical Casimir 
effect (DCE) also holds for a single mirror subject to 
nonuniform acceleration in empty space Although 
receiving considerable interest 0, Q since its theoretical 
discovery, there is still no experimental verification of the 
DCE. This is mainly due to the fact that the rate of pho- 
ton production is nonnegligible only when the mirror ve- 
locity approaches the speed of light, ruling out the use of 
massive mirrors [3| • Proposals for the experimental veri- 
fication of the DCE have suggested rapidly changing the 
field boundary conditions in other ways, e.g., using lasers 
to modulate the reflectivity of a thin semiconductor film 
[^, 01 or the resonance frequency of a superconducting 
stripline resonator [sj. 

Here we show that a coplanar waveguide (CPW) termi- 
nated by a superconducting quantum interference device 
(SQUID), as shown in Fig.[Tl is a very promising system 
for experimentally observing the DCE. The inductance of 
the SQUID can be controlled by the local magnetic flux 
threading the loop, giving a tunable boundary condition 
that is equivalent to that of a short-circuited transmission 
line with a tunable length. Because there is no massive 
mirror moving, the velocity of the effective boundary can 
approach the speed of light in the transmission line. The 
photon production from the vacuum can thus be made 
experimentally detectable, and its photon spectrum is 
identical to the one calculated in the DCE for a single 
oscillating mirror 

Building on work on superconducting circuits for quan- 




FIG. 1: (Color online) (a) A schematic diagram representing 
a CPW terminated to ground through a SQUID-loop. The 
effective inductance of the SQUID can be tuned by the mag- 
netic flux <I>ext(i), providing a tunable boundary condition, 
(b) The setup in (a) is equivalent to a short-circuited trans- 
mission line with a tunable length L^g, i.e., with a tunable 
"mirror" . We analyze this system using the input / output for- 
malism, which gives the spectrum of the scattered outgoing 
field, ^out, as a function of the incoming field, $in, in the 
coplanar waveguide. 



tum information [10|, [11|, [l^], there has recently been 
rapid progress in the field of circuit QED, where the in- 
teraction between artificial superconducting atoms and 
the electromagnetic field in microwave cavities is stud- 
ied. Recent achievements include strong coupling be- 
tween artificial atoms (qubits) and resonators [l3|, single- 
artificial-atom lasing [ij] , and Fock-state generation [l5| . 
Cavities with tunable frequencies and low dissipation 
have also been demonstrated [3, [I3|i ^^'^ has been 
shown that the resonance frequency can be changed by 
hundreds of linewidths on a timescale much shorter than 
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the photon hfetime in the cavity [18| . There are also re- 
cent theoretical suggestions to observe the DCE in cavity 
geometries including superconducting qubits 1^, 2^ , and 
also a suggestion to use a CPW where the center conduc- 
tor is replaced by an array of SQUIDs |2l| to simulate 
the Hawking radiation [23 |. 

We consider a superconducting CPW with character- 
istic capacitance Co and inductance Lq per unit length. 
The CPW is terminated at a; = through a SQUID-loop 
threaded by an external flux <i>oxt(0; shown in Fig.[T] 
Since the system contains Josephson junctions, it is con- 
venient to describe the electromagnetic field in the CPW 
line by its phase field ^{x,t) = J dt' E(x. t'), i.e., the 
time-integral of the electric field E{x,t) [23[. The phase 
field obeys the massless Klein-Gordon equation and, in 
second quantized form, is {x < 0): 



$(x,t) = 
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where the ai} operator annihilates a photon with 
frequency uj propagating to the right (left) with velocity 
V — \/ y/C^LQ and wave vector k^^ = lo/v, and satisfies 
the commutation relation [a™ (out)-jtj — _ 

Lo'). The characteristic impedance of the CPW is = 

We first consider a symmetric SQUID, where the two 
junctions have equal capacitances (Cj_i — Cj,2 — C/2) 
and Josephson energies {Ej^i = Ej^2 — Ej), and later 
comment on the effects of asymmetry. The SQUID effec- 
tively behaves as a single junction with a tunable energy 



Ejif) = Ej ^2 + 2 cosif), 



(2) 



where / = 27r<I>ext(i)/^o, and $o — h/2e. The junction 
can equivalently be characterized by its tunable (Joseph- 
son) inductance Lj{f) = ($o/27r)^/i?j(/), as long as 
the phase dynamics is slow compared to the plasma fre- 
quency Up if) = l/v/CLj(/), and the SQUID is only 
weakly excited. 

The effective boundary condition for the field imposed 
by the SQUID can be derived using quantum network 
theory . Starting from the classical Lagrangian for the 
circuit, the Heisenberg equations of motion are obtained 
from canonical quantization. For the system under con- 
sideration this analysis was performed in Ref . [25*1 . Here 
we are now interested in macroscopic SQUID junctions in 
the phase regime, i.e., when the charging energy is small 
compared to the Josephson energy, {2ef /2C ^ Ej{f), 
and the quantum fluctuations of the phase across the 
SQUID are small. In this regime, the boundary condi- 
tion at X = becomes 



idmt)^^d_^^,^ (3) 



where the last term can be neglected since we are consid- 
ering dynamics much slower than the plasma frequency 
of the SQUID. We note that this boundary condition de- 
pends parametrically on the tunable effective Josephson 
energy, Ej{f), providing the tunable boundary condition 
that will be essential for the remaining discussion. 

By inserting the field, Eq. ([T]), in the boundary condi- 
tion, Eq. ([3]), and performing an integration over time, 
we find the corresponding boundary condition in fre- 
quency space. For a static external flux /o, giving a static 
Josephson energy Ej = Ej{fo), we can solve for the out- 
put operators a°"* in terms of the input operators a™: 
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where the effective length 

-^cff = icff(/o) = ( ^ 
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is motivated by comparison to a short-circuited transmis- 
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sion line of length L and its phase factor —e^^^^. 
effective-length interpretation is valid for ki^L^f^ ^ 1, 
i.e., for frequencies where the SQUID effective length 
is smaller than the wavelength, or equivalently, lo <C 
ZoCujp, which is satisfied for the parameter regime that 
we are considering below. 

For a time-dependent external fiux, resulting in the 
Josephson energy Ej[f{t)] ^ Ej + 6Ej{t), we can write 
the solution in the form, 



J —oo 

X [e(c.)(aL", + + e(-c.)(aL", + a°V*)t] , (6) 



where, 
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and where O is the Heaviside step function. 

For a small-amplitude harmonic drive SEj(t) — 
SEj cos((jJdt), where SEj <C Ej, the effective length mod- 
ulation is also harmonic, Lcs{t) = L'^g + SL^g cos{uJdt), 
with amplitude SL^s = L'^si^Ej/Ej). In this case, we 
can evaluate the integrals in Eqs. ([B][7]) and perturba- 
tively solve the resulting equation for the output opera- 
tors in terms of the input operators, 
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VHu^-^d)\ [e(c. - c.,)aL"-., + e(c., - c.)(aL" -Jl) -(8) 

Considering a large-amplitude harmonic drive, we can 
expand Eqs. (07]) in terms of sideband contributions 

UJn = UJ + nUJd, 



dx 



N 

E 



c„ [eK)aL"„ + e(-a>„)«)^] , (9) 
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and numerically solve (for c„) the resulting set of linear 
equations to arbitrary order iV, where N is the number 
of sidebands included in the calculation. Note that we 
still require both the driving amplitude and the driving 
frequency to be small compared to the plasma frequency 
of the SQUID. Using this expansion we have numerically 
verified that a°"* and (a°"*)t obtained from Eq. (O sat- 
isfies the correct commutation relation [a°"*, (a™*)^] — 
S{lu — oj'), and, similarly, that the perturbative solution 
in Eq. ([8]) satisfy the same commutation relation to first 
order in the perturbation parameter 5Ej. 

Both the perturbative and numerical approaches give, 
in principle, all properties of the output field in terms of 
the input field, and they will be used below to calculate 
the output photon-flux density = ((a°"')ta°"*), as 
a function of mode-frequency oj, for the thermal input 
fields ~ l/(exp(?itj/fcsT) — 1). Using this photon-flux 
density, the number N of generated photons per second, 
in a bandwidth Atj, is given by 



N ^ — I duj 
2tt 



Au 



Alu 
2^ 



(10) 



where the approximation is valid for a small bandwidth 
Aw, where the relative change in is small. 
For a small-amplitude harmonic drive, we find 
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^(wd - t^) 0(^^d - t^), (11) 



where we have neglected terms containing the small fac- 
tor since we are considering fcsT <C ajd/27r. The 
output-field photon-flux distribution in Eq. (jlip can be 
decomposed into three components: The first two are 
of classical origin; elastically reflected thermal photons 
(first term) and up-converted thermal photons (second 
term). The third term is a purely quantum mechanical 
effect which originates from the vacuum fluctuations. We 
note that the spectrum of this quantum mechanical ra- 
diation is identical to the spectrum of the single-mirror 
dynamical Casimir effect 

The photon-flux-density spectrum of the quantum me- 
chanical radiation has a different frequency dependence 
compared to that of the reflected thermal photons and 
the two effects can therefore be clearly distinguished from 
each other. A signature of the quantum radiation is 
the parabolic shape in the photon-flux-density spectrum, 
which has a maximum at uJd/'^, whereas the photon-flux 
density for the reflected thermal fleld has maxima at zero 
frequency and at ujd, see Fig. [2] Furthermore, in this 
quantum mechanical radiation process, the photons are 
created in correlated pairs ((a| 



^ 0) with 

frequencies that sum up to the driving frequency, result- 
ing in a squeezing spectrum [26] with maximum squeez- 
ing at LOd/2. 
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FIG. 2: (Color online) The photon-flux density as a func- 
tion of mode- frequency uj. The red curve shows the classi- 
cal prediction, and the blue (analytical) and green (numeri- 
cal) curves show the quantum mechanical predictions. The 
numerical results were calculated using 10 sidebands below 
and above the center- frequency oj. The discrepancy between 
the analytical and numerical results is due to the higher- 
order sideband contributions that is only accounted for in 
the numerical calculations. Here we have used the parame- 
ters ujp/2-K ^ 36 GHz, = tJp/2, SEj = E^A, C = 90 fF, 
Zo ~ 55 f2, u ~ 1.2-10* m/s, and the input-field temperatures 
50 mK (solid lines) and zero Kelvin (dashed lines). 



Using the higher-order expansion in sideband contri- 
butions given in Eq. ([9]) we can write the photon-flux 
density as 



TV 



(12) 



n=-N 



where c„ are numerically-obtained coefficients. Each 
higher-order sideband gives an additional parabolic con- 
tribution, between zero and nujd, to the photon-flux- 
density spectrum. This explains the small discrepancy 
between the analytical and numerical results in Fig. [5) 

The experimental verification of the DCE in this sys- 
tem is possible, if for realistic circuit parameters the 
quantum contribution to the photon flux is distinguish- 
able from the classical thermal contribution. Figure [2] 
shows the photon-flux spectral density with the thermal 
and quantum contributions, for a moderate input mode 
temperature (50 mK) and ty pica l parameters for super- 
conducting electrical circuits [l8|. In this case, the DCE 
contribution to the photon flux is considerably larger 
than the classical contribution, in a wide frequency range 
between ^ LOd/'^ and LJd- The power per unit bandwidth 
in this frequency range gives an energy comparable to a 
few mK, which should be compared with the noise tem- 
perature of the typical amplifiers of a few K. Using lock-in 
techniques and long detection times, these power levels 
are clearly detectable. Plotting the radiation power per 
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FIG. 3: (Color online) The effective noise temperature, rgfj, 
in the output field of the CPW, versus the input-field temper- 
ature, T, for purely thermal radiation (red dashed) and for 
both thermal and quantum radiation (blue solid), at uJd/'2, 
half the driving frequency. For low enough temperatures, 
the quantum radiation is significantly larger than the thermal 
noise level. Here we used the same parameters as in Fig. |5] 



sentially identical in the classical and quantum cases. 

In conclusion, we have studied superconducting copla- 
nar waveguides with parametrically modulated boundary 
conditions and characterized the spectrum of the photons 
that are generated in this process. The system can be 
considered as a solid state analogue to quantum optical 
setups with moving mirrors, known to generate photons 
from vacuum fluctuations 2\. In the present setup, we 
show that a weak harmonic modulation of the boundary 
condition can result in broadband photon generation, i.e., 
the dynamical Casimir effect, and we estimate that it is 
feasible to detect this radiation in realistic experimental 
circuits. 

We would like to thank L. Tornberg, V. Shumeiko, S. 
Ashhab, G.J. Milburn and P. Bertet for helpful discus- 
sions. This work was in part supported by the NSA, LPS, 
ARO, and NSF Grant No. EIA-0130383. 



unit bandwidth (in temperature units) at half the drive 
frequency as a function of the input-field temperature 
(see Fig. ^ further illustrates how the quantum radia- 
tion dominates over the thermal radiation for sufficiently 
low temperatures. With these parameters (see the cap- 
tion of Fig. [2]) the crossover takes place around 70 mK, 
and the photon production rate A'' in a Aa;/27r = 100 
MHz bandwidth around 07^/2 is N ~ 10^ photons per 
second, see Eq. ([TU|) . 

We have so far only considered symmetric SQUID 
devices. Here we analyze the case where the SQUID 
junction capacitances (Cj^i, Cj^2) and Josephson energies 
{Eji,Ej2) are slightly asymmetric, i.e., Cj 3> ACj and 
Ef » AEj, where Cf = Cj,i+Cj,2, ACj = Cj^2-Cj,i, 
Ej = Ej^i + Ej^2, and AEj = Ej^2 ~ Ejs- Asym- 
metric capacitances give rise to a source term of the 
form iACj^/(t) in the boundary condition Eq. 
This source term produces a coherent oscillating response 
only at the driving frequency (neglecting the sidebands 
at zero and 2wd), which corresponds to a sharp addi- 
tional peak around ujd in the photon-flux-density spec- 
trum. The broadband feature below ujd is therefore unaf- 
fected by small asymmetries in the junction capacitances. 
An asymmetry in the Josephson energies gives rise to 
a correction to the effective SQUID Josephson energy, 
Ej{f), by the factor (1 - 2AEj/E^). For smaU asym- 
metries AEj <C E^, there will be a negligible reduction 
of the photon-flux density. 

Finally, we note that it would also be possible to study 
the DCE in a tunable cavity geometry J3|. However, 
here the task of clearly separating the DCE from the 
classical effect of parametric amplification is much more 
demanding since the stable states of the system are es- 
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